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Abstract 

We investigate a generalized Ising action containing nearest neighbour, next to nearest neighbour 
and plaquette terms that has been suggested as a potential string worldsheet discretization on cubic 
lattices by Savvidy and Wegner. We use both mean field techniques and Monte-Carlo simulations to 
sketch out the phase diagram. 

The Gonihedric (Savvidy- Wegner) model has a symmetry that allows any plane of spins to be 
flipped with zero energy cost, which gives a highly degenerate vacuum state. We choose boundary 
conditions in the simulations that eliminate this degeneracy and allow the definition of a simple fer- 
romagnetic order parameter. This in turn allows us to extract the magnetic critical exponents of the 
system. 
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1 Introduction 



In a series of recent papers [Q Savvidy, Wegner and co-workers suggested a Gonihedric random surface 
action which could be written as 
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5 = i E (1) 

<ij> 

on triangulated surfaces, where 6(aij) = |tt — a.ij\^ with C, some exponent and o;,j is the angle between 
the embedded neighbouring triangles with common link < ij >. This action is a robust discretization of 
the linear size of a surface, which is a well-defined geometrical notion that may be constructed in various 
equivalent ways. It was intended as an alternative to gaussian plus extrinsic curvature lattice actions of 
the form 

s = E(^-^) 2 + A E ( 2 ) 

<ij> A;,Aj 

which have been much explored Q as discretizations of rigid membranes and strings ||. Although a 
simulation showed that the action of equ.([I]) produced flat surfaces Q], potential problems arising from 
the failure to suppress the wanderings of vertices in the plane were pointed out in Q for the action with 
C = 1' Possible ways to cure this are to add additional Gaussian or linear terms to the action || or, more 
satisfactorily, to simply choose ( < 1. 

A study of the scaling of the string tension and mass gap in a dynamically triangulated model with 
an additional linear term produced inconclusive results [[7j, as have simulations of the gaussian plus 
extrinsic curvature action, because of the difficulties of simulating a complicated action on a dynamical 
surface. There is thus some incentive to investigate alternative approaches to such surface models where 
the computational costs are less onerous. One such approach for regularizing the Gonihedric action is to 
restrict the allowed surfaces to a (hyper)cubic lattice. This has been pursued in some detail analytically in 
H [| |l(| and one numerical simulation carried out in three dimensions |J . Recently Pietig and Wegner |jTTf 
have demonstrated with a Pcicrls contour argument that a transition does exist in the three-dimensional 
case and obtained a bound on the critical temperature that is not contradicted by the simulations. The 
crucial observation for this work is that the surface theory on a cubic or hypercubic lattice can be written 
equivalently as a generalized Ising action, where the boundaries between the spin clusters are the surfaces 
of the original model. The latticized Gonihedric model assigns a non-zero action only to right angled 
bends in the surface and self-intersections. Normalizing the weight for a right-angled bend on a link 
appropriately leaves a free parameter k that gives the relative weight of a self-intersection of the surface 
on a link. The energy of a surface on a cubic lattice is thus given explicitly by E = ri2 + Ann^, where ri2 
is the number of links where two plaquettes meet at a right angle and is the number of links where 
four plaquettes meet at right angles. 

A Hamiltonian which reproduces this energy on a cubic lattice has the form 

H = 2k E CTi °j ~ f E ° ia i + ~y~ E ( 3 ) 

<ij> [i,j,k,l] 

which is of generalized Ising form and contains nearest neighbour (< i,j >), next to nearest neighbour 
(<< i, j ») and round a plaquette k, I}) terms. Such actions are not new, having been investigated 
in some detail using both mean field methods and simulations in |T^ ]. However, the particular combination 
of coefficients arising in equ.(^) was not considered explicitly in this work because it corresponds to a 
particularly degenerate set of couplings. This degeneracy manifests itself in an extended symmetry in 
the model: for any value of n it is possible to flip a plane of spins with no energy penalty, providing 
the flipped plane does not intersect any existing portion of surface. This gives a vacuum degeneracy 
reminiscent of a gauge theory, the difference with a true gauge theory being that here the symmetries 
are quasi-global, giving only the freedom to flip entire planes rather than local spins. Related surface 
models have also been simulated directly in [fl3|| , but again the set of coefficients appearing in equ.(3) 
was not explored in this work. A very rich phase structure was observed in fl^| , in common with other 
Ising models with extended interactions []l4| of various sorts which display first and second order phase 
boundaries as well as incommensurate phases. Given the generic richness of the phase diagrams for such 
generalized Ising models and the additional symmetries present in the Gonihedric model, the action of 
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equ. (|3j) merits investigation from purely statistical mechanical considerations as well as from the point of 
view of finding potential continuum string theories. 

In the context of string theory one is looking for a continuous transition (or transitions) at which 
a sensible continuum surface theory may be defined. It is perhaps worth recalling that even this does 
not guarantee a good continuum surface theory. The interfaces in the standard nearest neighbour Ising 
model in three dimensions, which has a continuous phase transition, have been investigated in some detail 
recently and found to be very porous objects, decorated with lots of handles at the scale of the lattice 
cutoff jl5|. Ideally one might hope that the surfaces generated by the Gonihedric action were smoother, 
given that it is derived from a sort of stiffness term. 

Our motivation in this paper is to investigate the action of equ.(^|) in order to sketch out the gross 
features of its phase structure. We concentrate on k = 1 as in Q, but also discuss other values. A few 
cautionary words are in order before we go on to discuss the mean field approach and simulations. As 
we have noted the ground state of the action in equ. (J3J) is very degenerate as planes of spins parallel to 
any of the cube axes can be flipped at no energy cost. In the case n = the degeneracy is even larger 
as diagonal planes may be flipped now also. The ability to flip arbitrary spin planes makes defining a 
magnetic order parameter rather problematic. Even the staggered local order parameters defined in p2[ 
would miss the lamellar phases with arbitrary intersheet spacings that could be generated at no cost by 
flips of spin planes. 

The simulations for k = 1 in || were restricted to measuring the energy and specific heat as the 
exhaustive global order parameters suggested there 



M M = 




(where erf (vac) is one of the possible vacuum spin configurations with /i = 1,2...2 3L for k = or 
/i = 1, 2 ... 3 x 2 L for k ^ on a lattice of size L) would have been prohibitively slow to measure on even 
moderately sized lattices. It is possible to do rather better, however, by making use of the freedom in 
choosing boundary conditions on a finite lattice. It is customary to employ periodic boundary conditions 
in attempting to extract critical exponents from a simulation as these tend to minimize the finite size 
effects. It is clear that fixed boundary conditions in the Gonihedric model would penalize flipped spin 
planes by at least a perimeter energy, at the possible expense of greater finite size effects. A quick test 
simulation shows that such boundary conditions do pick out the purely ferromagnetic ground state from 
the many equivalent possibilities and allow the measurement of the simple ferromagnetic order parameter 




One can, in fact, have the best of both worlds by continuing to employ periodic boundary conditions to 
reduce the finite size effects whilst fixing any two perpendicular planes of spins to pick out the ferromag- 
netic ground state. The fixed planes of spin are, in effect, more akin to a gauge fixing of the spin flip 
symmetry than boundary conditions per se. In the simulations reported later in this paper we employed 
three fixed perpendicular planes of spins as a safety measure, with essentially identical results. 



2 Zero Temperature and Mean Field 

As the Gonihedric model is a special case of the general action considered in [[l2| we can apply the methods 
used there for both the zero temperature phase diagram and mean field theory. For the zero temperature 
case this involves writing the full lattice Hamiltonian as a sum over individual cube Hamiltonians 

K \ ^ k \ -\ 1 — K \ - 

h c = - °-i<rj }^ ai(J i + 1^ WjVkcn (6) 

<i,j> [i,j,k,l] 

and observing that if the lattice can be tiled by a cube configuration minimizing the individual h c then 
the ground state energy density is eo = min h c . 

We list the inequivalent spin configurations on a single cube and their multiplicities in Table. 1 using 
the same notation as [|l2| but with our choice of couplings to highlight the degeneracies that appear with 
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the Gonihedric action. In the list of spins the first column represents one face of the cube and the second 
the other. In the table two configurations are considered equivalent if one can be transformed into the 
other by reflections and rotations or if they are related by a global spin flip. The antiferromagnetic image 
of a configuration is obtained by flipping the three nearest neighbours and the spin at the other end 
of the cube diagonal from a given spin and is denoted by an overbar. With the Gonihedric values of 
the couplings the freedom to flip spin planes is clear even at this level as ipo, which would represent a 
ferromagnetic state when used to tile the lattice, and ipe which would represent flipped spin layers, have 
the same energy for any value of k. The higher energy configurations tp4 and iftg are also identical. The 
degeneracies increase when n = 0, the club of states of energy —3/2 is now composed of V'o > V'o > ^6 > V'g 
and various extra degeneracies appear for higher energy states. The new ground states ^oj^B reflect the 
additional freedom to flip diagonal planes of spins that is present at k = 0. 

In the mean field approximation the spins are in effect replaced by the average site magnetizations. 
The calculation of the mean field free energy is an elaboration of the method used above to investigate 
the ground states in which the energy is decomposed into a sum of individual cube terms. The next 
to nearest neighbour and plaquette interactions in the Gonihedric model give the total mean field free 
energy as the sum of elementary cube free energies (f>(m c ), given by 

<t>{mc) = - ^ m " lTa 3 + 4 mim i 

<i,j>dC <<i,j»CC 
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| ^2 mjmjmkmi + — y^[(l + mj)ln(l + rm) + (1 - mj)ln(l - m,)] (7) 

[iJ,k,l]CC iCC 

where mc is the set of the eight magnetizations of the elementary cube. This gives a set of eight mean- field 
equations 

Wpcd =Q (8) 

orrii ( i= i... 8 ) 

(one for each corner of the cube) rather than the familiar single equation for the standard nearest neigh- 
bour Ising action. More explicitly, we have 

mi = tanh[A(3K(m2 + + ms) — 2(3K(m^ + ttlq + mg) 
+2(3(1 — K)(m2m 3 m4 + mim^rn§ + m^rn^m^ 



ms = tanh[A(3K(m2 + + 1117) — 2/3«;(mi + m 3 + m%) 

+2/3(1 — K)(m,3m4m,7 + m\m^m,^ + msmemy)] (9) 

where we have labelled the magnetizations on a face of the cube counterclockwise 1 ... 4 and similarly 
for the opposing face 5 ... 8 as shown in Figure. I. If we solve these equations iteratively we arrive at 
zeroes for a paramagnetic phase or various combinations of ±1 for the magnetized phases on the eight 
cube vertices, and the mean field ground state is then give by gluing together the elementary cubes 
consistently to tile the complete lattice, in the manner of the ground state discussion. 

Turning loose a numerical solver on the mean field equs.(^|) gives generically a single transition to one of 
the phases listed in Table. I from the high temperature paramagnetic phase. The transition temperatures 
and the resulting low temperature phase are listed in Table. 2. We have taken the liberty of carrying out 
global flips where necessary to tidy up the table. Rather remarkably, we see that apart from k = the 
transition appears to be to the simple ferromagnetic phase, tJjq. However, remembering that V'o and tpe 
have the same energy the best we can say is that we end up in a layered phase with arbitrary interlayer 
spacing in all directions. Although the n = case appears to be superficially different, the ipQ phase that 
is found at low temperature here is one of the phases that is degenerate with ip^ and i/jq when k — 0. 
Although k = 1 fits the pattern as far as a transition to ^>o,6 at decreasing [3 is concerned it appears 
to be rather atypical in that further transitions are observed at larger /3. However, this is a numerical 
instability that is peculiar to this particular value of k. It was observed in jl2j that an iterative solution 
of the mean field equations written in the form 



in 



= f[E ti (m n )} (10) 
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where E is the individual cube Hamiltonian could fail to converge if an eigenvalue of dm\ n /dm™ was 
less than —1. Modifying the equations to 

m („ +1) = (f[E4m n )]+am?) 
1 1 + a 

for suitable a cures this. This is precisely what happens for k = 1, where introducing a non-zero a 
suppresses the extra "transitions" . 

In summary, the mean field theory suggests a rather simple phase diagram for the Gonihcdric model 
with action equ. (||) , with a single transition from a paramagnetic phase to a degenerate "layered" phase 
that is pushed down to ft = at large k. The low temperature phase is generically of the ^o,6 type, 
apart from k — where we see a tpQ g phase that is degenerate with these. The degeneracy of the ground 
states that are indicated by these results are, as they should be, consistent with the symmetries of the 
original full action. We now go on to see how the zero-temperature and mean field results tally with a 
direct Monte-Carlo simulation. 



3 Simulations 

We carried out simulations with n = 1 on lattices of size 10 3 , 12 3 , 15 3 , 18 3 , 20 3 and 25 3 and for k = 
0, 2, 5, 10 on lattices of size 10 3 , 15 3 , 20 3 and 25 3 . Unless stated otherwise periodic boundary conditions 
were imposed in the three directions and three internal perpendicular planes of spins fixed to be +1. We 
carried out 50K sweeps for most (3 values, increasing to 500K sweeps near the observed phase transition 
point. Measurements were carried out every sweep after allowing sufficient time for thermalization. 
A simple Metropolis update was used because of the difficulty in concocting a cluster algorithm for a 
Hamiltonian with such complicated interaction terms. The program was tested on the standard nearest 
neighbour Ising model and the some of the parameters used in the generalized Ising models of ||l2f to 
ensure it was working. 

We measured the usual thermodynamic quantities for the model: the energy E, specific heat C, 
(standard) magnetization M, susceptibility x ancl various cumulants. As the large (3 limit of the energy 
should be determined by the zero-temperature analysis of the preceding section, we consider the energy 
first. The absolute value of the energy for various n on lattices of size L = 20 is plotted against (3 in 
Figure. 2, where we can see that the zero temperature prediction of 3(1 + k)/2 is satisfied with good 
accuracy for sufficiently large (3. We can therefore observe that the zero-temperature/mean-field analysis 
has correctly identified the ground state(s) of the theory: -0 O ,6 for n ^ 0; or ^0,0,6,6 for k = as these are 
the only states with the observed energies. Having satisfied ourselves that the simulation is finding the 
correct ground state energy, we can go on to consider extracting some of the critical exponents for the 
transition. In what follows we will, as advertized, discuss in some detail the case n = 1 before commenting 
more briefly on the other values of k. 

With only periodic boundary conditions the possibility of the simple ferromagnetic ordered state i/'o 
can be excluded by looking at the magnetization M, which for all k is either zero or fluctuates wildly 
as p is changed, reflecting the freedom to flip spin planes. Curtailing this freedom by fixing the three 
perpendicular spin planes picks out the transition to a simple ferromagnetic ground state and the low 
temperature limit of the magnetization becomes one for all k. The magnetization cumulant 

< M 4 > , . 

Um = (12) 

is well defined once the three spin planes are fixed, but it does not show the standard behaviour of a low 
(3 limit of three and a high [3 limit of one, asymptoting to a value slightly larger than one at low [3, as can 
be seen in Figure. 3. This is because the fixed planes still leave a residual magnetization at low j3, which 
is sufficient to make this limit look magnetized for the sizes of lattice we simulate. Nonetheless, it is still 
possible to apply the usual scaling analysis in the critical region, and the crossing of the cumulant plots 
for different latttice sizes gives an estimate of f3 c = 0.44(1) for the critical temperature, which is in good 
agreement with the value reported in [^) that was extracted by looking at the change in behaviour of the 
spin/spin correlation as (3 C was approached. As noted in || this [3 C is very close (in our case within the 
error bars) to that of the standard two-dimensional Ising model on a square lattice. 
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The rather sharp nature of the crossing, or more accurately collapse down to a single line, makes it 
difficult to extract a value for v from the ratio of the slopes of the Binder's cumulant curves at the critical 
point, so we take a different tack and consider the scaling of the maximum slope, which we would also 
expect to scale as L 1 /" []. This gives an estimate of v = 1.2(1). We can also look at both the finite size 
scaling and direct fits to the susceptibility x, namely x — AL» and \ — A\j3 — /3 C | 7 , to attempt to extract 
v. We choose this in preference to the specific heat fits C = B + DL^ and C = B + D\(3 — (3 c \ a because 
of the absence of an adjustable constant. The susceptibility data is plotted in Figure. 4, and shows a clear 
peak. This should be contrasted with the case of no fixed planes where the (3 > (3 C region is rendered 
meaningless noise by the lack of a well-defined magnetization for the myriad of ground states. We find 
the finite size scaling fit gives 7/1/ = 1.79(4) with a high quality, and feeding the critical value of (3 C = 0.44 
into the direct fit on the L = 25 lattice gives 7 = 1.60(2) with rather lower quality. The deduced value 
for v is thus 1.10(5). These fits give values close to those for the standard two-dimensional Ising model 
with only nearest neighbour interactions, where we have 7 = 1.75, v = 1. 

As a consistency check on our values of (3 C and v, we plot the (3 values where the specific heat peaks 
and the f3 values where the maximum slope of the Binder's cumulant curve occurs, both of which can 
serve as estimates of the pseudocritical temperature on finite size lattices, against L~~ . We would expect 
this to be a straight line with intercept (3 C . The plot is shown in Figure. 5 for the choice v — 1 with 
other values in this neighbourhood giving essentially identical results. We can see that the estimate of 
j3 c = 0.437(7) coming from the two possibilities is both self-consistent and in agreement with the value 
obtained from cumulant crossing. 

The above, apparently consistent, set of results presents us with something of a dilemma when it 
comes to the analysis of the specific heat peak, which is shown in Figure. 6. The hyperscaling relation 
a = 2 — vd indicates that, if a value of v ~ 1 is to be believed, the specific heat should display a cusp 
(a = 2 — vd ~ — 1) rather than a divergence. This does not appear to be the case for the data in the 
figure. Setting aside the hyperscaling result for the moment and performing direct power law fits to 
C = B + DL^ and C = B + D\(3 - f3 c \ a gives poor fits. A logarithmic fit of the form C = B + D \og{L) 
or C = B + D log(/3 — (3 C ) gives much better, but still not particularly good, results so the evidence is 
inconclusive. 

Another line of attack for obtaining an estimate of a is to use the finite size scaling of the energy 
itself E ~ Eq + E\L~~z~ . With direct measurements in this form one gains nothing in general over the 
specific heat fits as there is still a regular term Eq to be dealt with, although for models with higher than 
second order transitions this approach may be preferable Jig] , However, if one has measurements for two 
different sets of boundary conditions available the regular term would be expected to be the same for both 
and the energy difference could be used for a simple power law fit to extract (a — We are currently 

measuring the string tension in the Savvidy model using two sets of fixed boundary conditions p7| , one 
with all positive spins and one with half positive and half negative spins 0. For these measurements we 
would expect 

AE = E ++ - E+- = AL^ (13) 

where E ++ is the energy for all positive spins and E^ is the energy for half positive and half negative 

spins. A fit gives (a — \)/v — —1.3(2), which is still marginally consistent with a = 0. 

There are two possible interpretations of the results. The first is that the value of v measured is 
simply wrong, not inconceivable as it appears either as a derived quantity from the slope of the cumulant 
or from the two fits to the susceptibility rather than being measured directly. However, a second possible 
interpretation is to accept the fits to v at face value and posit that the model sees an effective dimension 
of d = 2 in order to avoid violating the hyperscaling relation. In this case we could recover the full 
set of two-dimensional Ising model exponents. Although this would be highly unusual, it should be 
remembered that the energy in the model is essentially linear in form rather than being an area, so there 
is some resemblance to a two-dimensional model. 

A direct fit to the magnetization exponent on the largest lattice size M ~ \(3 — (3^ (with apologies 
for the profusion of betas!) with j3 c fixed to be 0.44 gives j3 — 0.12(1), but the quality is low, whereas a 
finite size scaling fit gives a much lower value of f3/v = 0.04(1). It is possible that the fixed spin planes, 
whose residual magnetization we have not allowed for in the fits, are biasing the finite size scaling fit, but 

1 This tactic works well in, for instance, simulations of Ising models coupled to two-dimensional quantum gravity. 
2 It is not possible to use antiperiodic boundary conditions in the Savvidy model to enforce an interface because of the 
plane flip symmetry. Something more coercive, in the form of these fixed boundary conditions is required. 
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intuitively one would expect their effects (of order 3/L) to increase rather than decrease the estimated 
exponent by pushing up the measured magnetization on the smaller lattices. 

We now discuss more briefly the other k values that were simulated, namely k = 0, 2, 5, 10. Firstly 
we can note that in qualitative terms the transition appears similar to the k = 1 case, a not entirely 
trivial result as /t ^ 1 introduces round a plaquette interactions that are missing for k = 1. With the 
periodic plus fixed plane boundary conditions we still have peaks in the susceptibility and specific heat 
and a ferromagnetically ordered phase appearing at low temperature. The mean field analysis suggests 
that as k is increased (3 C should drop sharply. This is not observed in the simulations, the crossing of the 
Binder's cumulants indicating no change within the error bars for the estimates of C , giving (3 C — 0.44(1) 
from k = 1 to k = 10. There is a sharper difference with the k — results which show a crossing at 
f3 c = 0.50(1). In general mean field theory will underestimate /3 C , so the measured results are in agreement 
with this and do not contradict the bound obtained in [[llj . 

The similarity of the transitions for different k extends beyond (3 C . The measurements of j/v listed 
in Table. 3 below for all the non-zero k suggest that the critical behaviour is unchanged by varying k. 



K 


1 


2 


5 


10 


l/v 


1.79(4) 


1.6(1) 


1.9(1) 


1.75(6) 



Table 3: Fits to 7/f for the non-zero K values. 
The specific heat curves and magnetization present a similar story, with all looking similar to the k = 1 
case. From this evidence it would seem that varying k, at least for k > 1, gives little if any change in the 
exponents and transition temperature. 

The story is slightly different for k — 0. As we have already noted f3 c = 0.50(1), and other differences 
are apparent. Without the fixed spin planes (ie with only periodic boundary conditions) the susceptibility 
when P < P c for non-zero k values is similar to the fixed plane case described above and becomes 
meaningless when /3 > (3 C where the magnetization is ill-defined. The n — model presents qualitatively 
different behaviour in that the susceptibility is one for (3 < /3 C (— 0.5) and zero for (3 > (3 C . However, this 
step function behaviour disappears when the fixed spin planes are introduced and we recover a divergent 
peak as for the other k values. The phase structure for k = with the fixed spin planes also appears to 
be broadly similar to other n, giving a single transition to a low temperature magnetized phase. It would 
be interesting to examine in detail values of n between zero and one to see if there was a smooth change 
in, for example, f3 c as k — > 0. This would give some indication of whether k = really was a special 
point, or joined on smoothly to the continuum of non-zero values. A test simulation at k = 0.5 still gives 
very similar results to k — 1, for example. 

4 Conclusions 

We have conducted zero-temperature, mean-field and Monte-Carlo investigations of the generalized Ising 
model action suggested in || ^|, as a cubic lattice discretization of the Gonihedric string action Q 
using essentially the methods of |12f| . Although the phase structure of such generalized Ising models is 
generically very rich |]l4|| , the one parameter family of models examined here seems to be a fairly simple 
"slice" of the phase diagram, with one transition to a layered ground state when periodic boundary 
conditions are imposed. This degenerate layered state is a consequence of the plane spin flip symmetry 
that is present in the model for all k, but a judicious choice of boundary conditions - fixing enough 
perpendicular spin planes - allowed us to pick out an equivalent ferromagnetic ground state for the 
purposes of simulations. The zero-temperature/mean-field analyses are in agreement with the Monte- 
Carlo simulations on the nature of the ground state and its energy, but the simulations indicate a transition 
temperature that changes little, if at all, from its value at n = 1 (f3 c ~ 0.44) for other non-zero k values. 
The mean field theory on the other hand gives a fairly sharp decline in f3 c as k is increased. 

The simulations at k = 1 indicate that the fitted exponents, with the exception of the finite size scaling 
fit to /3/v, and even the critical temperature are all in the vicinity of those for the two-dimensional Ising 
model, though given our modest statistics it would be foolhardy to claim they were identical on the 
basis of the current fits. Comparison with the other non-zero n values also gives similar exponents and 
critical temperatures. There is some evidence that the n — model is a special case: - in the zero 
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temperature and mean field analyses more ground states are allowed and in the simulations a different 
critical temperature is observed and the behaviour of the susceptibility is radically different when no spin 
planes are fixed. 

An immediate extension of the current work, given the closeness of the fitted exponents to the two- 
dimensional fsing model, is to carry out a higher statistics simulation near the transition point in order 
to pin down the various exponents and (3 C more accurately. A further test of the critical behaviour of 
the model would be to investigate the scaling of the string tension as one approached the critical point, 
along the lines of ]17| . The various higher dimensional generalizations that were formulated in § |, |§ 
also merit investigation. 

If we return to our original stringy motivation it would be useful to characterize the surfaces generated 
by the Gonihedric action in the style of Jl5| to see whether they were any less "spongy" than those in the 
standard 3D Ising model. As a playground for exploring plaquette discretizations of string and gravity 
inspired models the generalized Ising models clearly have some interesting quirks that are worthy of 
further exploration. It would certainly be amusing to show that a candidate discretized string model was 
a close relation of the too-dimensional Ising model. 



5 Acknowledgements 

R.P.K.C. Malmini was supported by Commonwealth Scholarship SR0014. 



7 



References 



[1] R.V. Ambartzumian, G.S. Sukiasian, G. K. Savvidy and K.G. Savvidy, Phys. Lett. B275 (1992) 99; 
G. K. Savvidy and K.G. Savvidy, Int. J. Mod. Phys. A8 (1993) 3393; 

G. K. Savvidy and K.G. Savvidy, Mod. Phys. Lett. A8 (1993) 2963. 

[2] S. Catterall, Phys. Lett. B220 (1989) 207; 

C. F. Baillie, D. A. Johnston and R. D. Williams, Nucl. Phys. B335 (1990) 469; 
R. Renken and J. Kogut, Nucl. Phys. B354 (1991) 328; 
S. Catterall, Phys. Lett. B243 (1990) 121; 

C. F. Baillie, R. D. Williams, S. M. Catterall and D. A. Johnston, Nucl. Phys. B348 (1991) 543; 
S. Catterall, D. Eisenstein, J. Kogut and R. Renken, Nucl. Phys. B366 (1991) 647; 
J. Ambjorn, J. Jurkiewicz, S. Varsted, A. Irback and B. Petersson, Phys. Lett. B275 (1992) 295; 
J. Ambjorn, J. Jurkiewicz, S. Varsted, A. Irback and B. Petersson, Nucl. Phys. B393 (1993) 571; 
M. Bowick, P. Coddington, L. Han, G. Harris and E. Marinari, Nucl. Phys. B394 (1993) 791. 

[3] A. Polyakov, Nucl. Phys. B268 (1986) 406; 

H. Kleinert, Phys. Lett. B174 (1986) 335; 
W. Helfrich, J. Phys. 46 (1985) 1263. 

[4] C.F. Baillie and D. A. Johnston, Phys. Rev. D45 (1992) 3326. 

[5] B. Durhuus and T. Jonsson, Phys. Lett. B297 (1992) 271. 

[6] C.F. Baillie, D. Espriu and D. Johnston, Phys. Lett. B305 (1993) 109. 

[7] C.F. Baillie, A. Irback and W. Janke and D. A. Johnston, Phys. Lett. B325 (1994) 45. 

[8] G. K. Savvidy and F.J. Wegner, Nucl. Phys. B413 (1994) 605; 
G. K. Savvidy and K.G. Savvidy, Phys. Lett. B324 (1994) 72; 

G. K. Savvidy, K.G. Savvidy and P.G. Savvidy, "Dual Statistical Systems and Geometric String," 
frep-th/9409031| ; 



G. K. Savvidy and K.G. Savvidy, Phys. Lett. B337 (1994) 333; 
G. K. Savvidy, K.G. Savvidy and F.J. Wegner, Nucl. Phys. B443 (1995) 565. 

[9] G. K. Bathas, K. G. Floratos, G. K. Savvidy and K.G. Savvidy, "Two Dimensional and Three 
Dimensional Spin Systems with Gonihedric Action" , hcp-th/9504054 . 

[10] G. K. Savvidy and K.G. Savvidy, "Interaction Hierarchy: Gonihedric String and Quantum 
Gravity," [hcp-th/9506184 . 

[11] R. Pietig and F. Wegner, "Phase Transition in Lattice Surface System with Gonihedric Action", 
Heidelberg preprint, December 95. 

[12] A. Cappi, P. Colangelo, G. Gonclla and A. Maritan, Nucl. Phys. B370 (1992) 659. 

[13] T. Sterling and J. Greensite, Phys. Lett. B121 (1983) 345; 
M. Karowski and H. Thun, Phys. Rev. Lett. 54 (1985) 2556; 
M. Karowski, J. Phys. A19 (1986) 3375. 

[14] W. Selke, Phys. Rep. 170 (1988) 213; 

D.P. Landau and K. Binder, Phys. Rev. B31 (1985) 5946. 

[15] V. Dotsenko, M. Picco, P. Windcy, G. Ha rris, E. Martinec and E. Marinari, "Self- Avoiding Surfaces 
in the 3-D Ising model" , |hcp-th/9504"07^ ; 

V. Dotsenko, M. Picco, P. Windey, G. Harris, E. Martinec and E. Marinari, Phys. Rev. Lett. 71 
(1993) 811. 

[16] W. Janke and C. Holm, J. Phys. A27 (1994) 2553. 

[17] M. Baig, D. Espriu, D. Johnston and R.P.K.C. Malmini, "String Tension in Gonihedric Ising Models" , 
to appear. 



8 



TABLES 
Table. 1 



State 


Top 


Bottom 


Energy 


Multiplicity 


ipo 


+ + 
+ + 


+ + 
+ + 


-3/2-3k/2 


2 


i>o 


- + 
+ - 


+ - 
- + 


-3/2 + 21k/2 


2 


in 


+ + 
- + 


+ + 
+ + 


-3k/2 


16 


V'i 


- + 
+ - 


+ + 
- + 


9k/2 


16 


^2,2 


- + 
+ - 


+ + 
+ + 


1/2 + k/2 


24 


^3 


+ + 


+ + 
+ + 


-1/2 - 3k/2 


24 


ip2 


- + 


+ + 
+ - 


-1/2 + 5k/2 


24 


^4 


- + 
+ + 


+ + 

+ - 


3/2-3k/2 


8 


Ipl 


- + 


- + 

+ + 


3/2-3k/2 


8 


^5 


- + 


+ + 
+ + 


-3k/2 


48 




- + 
+ 




k/2 


48 






+ + 
+ + 


-3/2 - 3k/2 


6 




■ + 


+ + 


-3/2 + 5k/2 


6 


VV,7 


- + 


+ - 
+ + 


1/2 - 3k/2 


24 



Table. 1: The inequivalent spin configurations of a single cube and the associated energies and degeneracies 
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Table.2 





13 
Pc 






0.0 


0.325 


+ - 
- + 


- + 
+ - 


0.25 


0.31 


+ + 
+ + 


+ + 
+ + 


0.5 


0.278 


+ + 
+ + 


+ + 
+ + 


1.0 


0.167 


+ + 
+ + 


+ + 
+ + 


2.0 


0.09 


+ + 
+ + 


+ + 
+ + 


5.0 


0.0335 


+ + 
+ + 


+ + 
+ + 


10.0 


0.02 


+ + 
+ + 


+ + 
+ + 


15.0 


<0.02 


+ + 
+ + 


+ + 
+ + 



Table.2: The ground state configurations and transition temperatures for various k. 
The states shown appear above the quoted temperature. 
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Figure 1: The labelling of the cube vertices for the mean field equations. 
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Figure 2: The energies for various k, all on lattices of size L = 20. 
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Figure 3: Binder's cumulant for n = 1 
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Figure 4: The susceptibility for k = 1 
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Figure 5: The pseudocritical temperature estimated from the specific heat and Binder's cumulant vs L 
at k = 1 
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Figure 6: The specific heat for n = 1. 
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